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In superconductors, and in other systems with a local U(l) gauge invariance, 
there are two mechanisms that form topological defects in phase transitions. 
In addition to the standard Kibble mechanism, thermal fluctuations of the 
magnetic field also lead to defect formation. This mechanism is specific to 
local gauge theories, predicts a qualitatively different spatial defect distribu- 
tion and is the dominant source for topological defects in slow transitions. I 
review the arguments that lead to these conclusions and discuss possibilities 
of testing the scenario in superconductor experiments. 
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1. INTRODUCTION 

The Kibble mechanism! gives a very simple and elegant description for 
the formation of topological defects in symmetry breaking phase transitions, 
and has been studied intensively in experiments with liquid crystals! and 
liquid Helium£K§El Even though the results do not agree completely with the 
theoretical predictions, the mechanism itself is still believed to be correct. 
Instead, the discrepancies are presumably due to the extra assumptions one 
has to make about the dynamical time and length scales of the system in 
order to derive quantitative predictions H! 

More recently, defect formation has also been studied in superconduc- 
tors! and Josephson junctions.! These systems differ from liquid Helium in 
a significant way: the order parameter, i.e., the Cooper pair, is electrically 
charged, and the symmetry that is broken in the phase transition is the local 
U(l) gauge invariance associated with electromagnetism. This means that 
the arguments behind the Kibble mechanism cannot be applied directly, and 
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there has been some discussion over the years on how it should be modified.0 
However, it was recently pointed out by myself and Hindmarshllil that 
it is not enough to simply modify the Kibble mechanism, because there is 
another, totally separate effect that leads to defect formation. While the 
Kibble mechanism is based on the dynamics of the order parameter field, 
this new scenario relies on thermal fluctuations of the magnetic field. It also 
leads to predictions that are qualitatively different from those of the Kibble 
scenario. 

Unfortunately, the experimental setup in Ref. |8] was such that the con- 
tribution from our mechanism was negligible. Basically, one would have to 
be able to detect the positions and signs of individual vortices in order to be 
able to test the mechanism. 



The aim of this paper is to review the mechanism of Ref. 11 and to 
discuss its predictions for a superconducting film. 



2. KIBBLE MECHANISM 

Let us first review the Kibble mechanism and consider a simple two- 
dimensional toy model with an U(l) symmetry, i.e., our order parameter ijj 
is a complex scalar field. The thermodynamics of the model is described by 
the free energy,H3 which we approximate in the standard way by 



d 2 x 



Vip* ■ S/ip + mVV + A (ip*ifj) 2 , (1) 



where m 2 ~ T — T c . 

When the critical temperature is approached from the symmetric phase, 
the correlation length of the order parameter, i.e., the coherence length £, 
diverges but it cannot grow arbitrarily fast, because at least it is constrained 
by the speed of light. Therefore, if the phase transition takes place in a finite 
time, the coherence length reaches a finite maximal value £, and if we divide 
the system in domains of radius £, the phase angle of the order parameter 
is uncorrelated between different domains. At points where three or more 
domains meet, there is a fixed probability that the phase angle cannot be 
smoothly interpolated between the domains, in which case a vortex is formed. 
Thus, the number density of vortices must behave as 

"Kibble ~ i~ 2 - (2) 



Because the mechanism is based on very general assumptions only, this 
result applies to all systems in which a global U(l) symmetry breaks down 
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spontaneously, for instance to 4 He. It can also be generalized to other sym- 
metry groups very easily. However, in order to make contact with experi- 
ments, one would have to be able to predict the value of £. An upper bound 
can be obtained from causalityjll but a more precise estimate requires extra 
assumptions about the dynamics of the system.i'i 



3. GAUGE INVARIANCE 

If the order parameter field is electrically charged, as it is, e.g., in the 
case of superconductors, it couples to the electromagnetic vector potential. 
Then the analogue of the approximate free energy ([fl) ist3 



F[ip,A] = J d 2 x 



^B 2 + Dip* ■ Dip + m 2 ip*ip + A (ip*ipy 



(3) 



where B is the magnetic induction,[]i? = d x A y —d y A x , and D is the covariant 
derivative D = V+ieA, and A is the electromagnetic vector potential (gauge 
field). 

The coupling to A has the important consequence that the U(l) sym- 
metry becomes local; the free energy is invariant under all transformations 

ip(x) -> exp(iA(x))ip(x), A(x) -> A(x) - -VA(f), (4) 

e 

where A(x) can be any real function of x. 

In particular, this local gauge invariance means that the phase angle of tp 
does not have any physical meaning, because it can always be transformed 
away by a suitable choice of A. For this reason, one cannot talk about 
domains of correlated ip anymore. In the global case, the gradient term in 
the free energy (|I|) tends to make the order parameter aligned, but in the 
local case, it is the covariant derivative Dip and not the gradient that will 
be minimized. In the special case in which A = 0, there is no difference, 
but in practice there are always non-zero thermal fluctuations present in 
the A field. If A is non-zero, the covariant derivatives cannot typically be 
minimized everywhere at the same time, and this leads to vortex formation. 

In fact, this effect can easily be seen in simple cases in which the initial 
value of the vector potential A has been set by hand. As a very simple 
example, let us think of a phase transition in the presence of an homogeneous 
magnetic field B(x) = Bq. When the system enters the broken phase, the 
dynamics tries to minimize Dip. Let us write ip(x) = v(x) exp(i6(x)), and 
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consider a circular curve of radius R. At every point on this curve, the system 
tries to minimize the covariant derivative, but that leads to the result 

/ dx,V9 = - I dx-A= - I d?xB= 2 -^. (5) 
Jc e Jc e Js e 

This means that if RBq > e, non-zero winding number is preferred, and con- 
sequently a vortex is formed. This formation of an Abrikosov vortex lattice 
is a well-known phenomenon in superconductor physics. 



4. PLANE WAVE 

In order to understand the dynamics of phase transitions in which the 
magnetic field is inhomogeneous, it is instructive to consider first a plane 
wave. Let us assume that there is a standing electromagnetic wave trapped 
inside our two-dimensional system, 

B(x,y) = -E?osin(/cx), (6) 

where k is the wave number of the wave. 

The free energy @ by itself determines only the thermodynamics of the 
model, not the dynamics. For that, we need the equations of motion, and 
here we will consider a simple, "relativistic" model, 

doEi = eijdjB + 2eImV>* AV' > 

V-E = 2elmip*d i/j, (7) 

where V(ip) = m 2 \tp\ 2 + A|^| 4 . Note that there is no external magnetic field 
acting on the system; instead, Eq. (|6|) merely specifies the initial conditions 
for the time evolution. 

Let us now "quench" the system into the broken phase by changing m 2 
rapidly. To be specific, let us assume that 

m 2 (t) = -ml—. (8) 

If tq is small enough in comparison to the the wave number k, the situation 
in the regions in which B is, say, positive looks exactly the same as in the case 
of a homogeneous magnetic field. Therefore, vortices with a positive sign will 
form. Likewise, in the regions with negative B, vortices with a negative sign 
form. Fig. [l] shows the configuration of the B field at different times during 
such a transition. In the initial state, the order parameter field ip was given 
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Fig. 1. The time evolution of the local magnetic induction B(x) in a 
slow phase transition, starting from a plane wave configuration. Shades 
from white to black indicate field values from positive to negative. Initially, 
the amplitude of the field is very small, but the the contrast of the plot is 
adjusted accordingly to make it visible. In the final state, the flux is confined 
into vortices inside which the value of B{x) is very high. The signs of the 
vortices are strongly correlated with the initial value of B at each point. 

tiny random fluctuations, and then the mass parameter was decreased so 
that the system underwent a phase transition to the broken phase. The plot 
shows clearly how the signs of the vortices are correlated with the initial 
value of the magnetic induction at each point. Note, however, that the 
agreement is not perfect. These vortices with "wrong" signs are formed by 
the Kibble mechanism. 

On the other hand, if tq is large enough, the transition is adiabatic. The 
amplitude of B decreases, and the system ends up in the equilibrium state, 
which contains no vortices because of the absence of any external magnetic 
fields. This shows that for each k there is a critical value of tq so that if the 
transition is faster, vortices will form. We can also invert this relation and 
define for each tq a critical wavenumber /c(tq). 



5. THERMAL EQUILIBRIUM 

In practice, we are not interested in creating vortices by external mag- 
netic fields or by standing electromagnetic waves. Instead, we would like 
to understand what happens when the transition starts from thermal equi- 
librium. However, similar arguments apply to that well, because 
thermal fluctuations alone can act as seeds for vortex formation. 

In thermal equilibrium, the state of the system is not described by a 
single field configuration, but by an ensemble of configurations in which the 
probability of any given configuration is proportional to exp(—E[ip,A]/T). 
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Here E[ip, A] is the energy of the field configuration, and is essentially equal 
to the free energy @. Because of the smallness of the gauge coupling con- 
stant e, we can neglect the coupling to ip and approximate the energy by 

E[A] » J d 2 x l -B{xf. (9) 

Thus the probability distribution of B is almost Gaussian, and the field 
strength is totally uncorrelated between all the points in the coordinate 
space. If we write the energy in the Fourier space, we find that the same is 
true there as well, 

i.e., the amplitude of each Fourier mode independent of all the others and has 
a Gaussian distribution with a fe-independent width. As long as the system is 
in the symmetric phase, each mode also behaves dynamically independently 
of the other modes. Therefore we can estimate that in a transition with 
quench time tq, the modes with |A;| < A;(tq) equilibrate, but the modes with 
\k\ > k{TQ) fall out of equilibrium and their distribution remains the same 
as it was in the symmetric phase before the transition. 

To make comparison with the Kibble mechanism easier, we denote the 
wavelength corresponding to k by £ = 2ir/k. The above implies that if we 
concentrate on a region S of radius £, the total magnetic flux through this 
region does not change significantly in the transition. Inside the region, the 
flux rearranges itself in order to minimize the energy, i.e., it forms a vortex 
"lattice" in which all the vortices have equal sign. Therefore, we can estimate 
the number of vortices in the region S after the transition by calculating the 
typical flux through it at the time of the transition, and this is approximately 
the same as it is in the symmetric phase. Because of the Gaussianity of the 
probability distribution, it is easy to calculate, 

(V/ p ) 2 » ($|) = J d 2 xd 2 x'{B{x)B{x')) = J d 2 xT = iri 2 T. (11) 

Therefore, the typical number of vortices through the region S in the broken 
phase is 

N S = 1>f p /*o « ^T 1/2 L (12) 
and, consequently, the number density of vortices is 

n « § « f T^r*. (13) 



Note that this estimate was for a two-dimensional system, and therefore 
e is not a dimensionless number like it is in three dimensions, but instead, 



Formation of Topological Defects 




Fig. 2. The winding number of a circular curve of radius r centred at a 
positive vortex, No(r), behaves differently in the Kibble mechanism (left) 
and in the local gauge theory (right). 

it has a dimension of energy 1 / 2 . The three-dimensional case has also some 
other differences, and it will be discussed in Section [F] 

The essential point in Eq. ( |l3|) is the relation n oc which differs 
from the prediction of the Kibble mechanism n oc £ -2 in Eq. (|2|). In order 
to really predict the vortex number density in either case, we must be able 
to calculate £ and that depends sensitively on the dynamics of the system. 
However, if we assume that the critical wavelength £ is or the same order 
of magnitude in both cases, Eq. (|l3l ) becomes greater than Eq. (||) when 
£> 2ir/eT 1 l 2 . Because long wavelengths correspond to slow transitions, we 
can conclude that this mechanism dominates in slow transitions. 



6. SPATIAL DISTRIBUTION 

Even without any assumption on £, our mechanism can be tested by 
studying the spatial distribution of vortices. One particularly simple way 
of seeing this is by considering the average winding number along a circular 
curve of radius r centred at a vortex of positive sign (see Fig. ^). In practice, 
we pick a vortex, draw a circle of radius r around it and count the number 
of vortices of the same sign as the one at the centre minus the number of 
vortices of opposite sign, and average this over all vortices. We denote this 
quantity by Nc{r). Obviously, iVc'(O) = 1, because then only the vortex 
itself is counted. 

Let us first consider the Kibble mechanism. In this case, the vortex 
number is determined by the phase angle 7 = arg ip of the order parameter 
field, and we can write 

N c (r) = { dr- V7. (14) 
Jc 
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The vortices are formed at the points where domains meet, and therefore 
the distance to the nearest vortex is roughly the domain size £. Therefore, 
as long as r <C £, Nc(r) does not change significantly as we increase r. 

However, when every point on the curve C is uncorrelated with 

the centre of the curve, because the correlation length £ is less than r. There- 
fore the integral (|l4|) must be independent of whether there is a positive or 
a negative vortex at the centre. This implies that it cannot have a preferred 
sign, and therefore it must vanish on average. Summarizing, we find that in 
the Kibble mechanism 



Nc(r)*{ l ' '~V (15) 




In our mechanism, the vortex number is determined by the magnetic 
field, and therefore it is more convenient to writd^ 

N c (r) = t d 2 xB(x), (16) 
JS(C) 

where S(C) is the region bounded by the curve C. In this case, the vortices 
are inside domains and all the vortices in a given domain have equal sign. 
When we start increasing r from zero, each new vortex we pick up has the 
same sign as the one at the origin, as long as r <C £. This implies that 
Nc(r) increases with increasing r. However, when r 3> £, the sign of any 
new vortex encountered is independent of the vortex at the centre, and 
therefore, on average, Nc{r) remains constant. Thus we conclude that 

N c (r)*{ 1 + Cr2 ' , . r ~f (17) 
I constant (> 1), r>£, 

where c is some constant. 

The two predictions (|l^) and ( |l~7|) are qualitatively different, and there- 
fore it is easy to determine which of the mechanisms dominates if one is 
able to measure the spatial distribution of vortices in the final state. In 
numerical simulations, this is straightforward, and we have carried out this 
measurement .tH 

In our simulations, we used a thin (120x120x5) three-dimensional lat- 
tice, and solved numerically the the equations of motion 

V-E = 2elmi/)*d ip. (18) 



In a suitably chosen gauge, Eq. 



arguments that lead to the prediction ( L5I) are not 



is valid in the gauge theory as well, but the 
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Fig. 3. The average of Nc{r) as a function of r after slow (tq 
fast (tq 



500) and 



10) transitions (from Ref. 11). The solid line is a benchmark 
curve, which shows the behaviour in the case where there are no correlations 
at all between vortices. In both cases the data points are significantly above 
this curve at short distances, which indicates positive correlations at short 
distances, just like predicted by Eq. (0). The signal gets stronger in slow 
transitions. 



The time step was 5t = 0.05 in units where the lattice spacing is one. The 
coupling constants were e = 0.3 and A = 0.18. Because A > e 2 , this cor- 
responds to a type-II superconductor. We first thermalized the system to 
temperature T = 6 in the symmetric phase using a hybrid Monte Carlo al- 
gorithm. Then, we evolved the system according to the equations of motion 
(|j~8|), changing the mass term as 

m 2 (t) = m 2 , — 5m 2 \ — arctan 1 — | . (19) 

V 37r T Q 3 / 

Here tUq = 1.4 is the initial value at t = — tq when the quench begins and 
5m 2 = 3.2 was chosen in such a way that the transition takes place at t ~ 0. 
In the final state, we measured the function Nc(r) and averaged it over 
several initial configurations taken from the thermal ensemble. The results 
are shown in Fig. || for two different quench rates. The solid line in the plot 
is to illustrate the effects of the finite system size. Because the total flux 
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through the system is fixed to zero by the boundary conditions, Nc(r) would 
not be unity even for a random distribution of vortices. Instead, it would 
decrease as Nc(r) = 1 — nr 2 /A, where A = 120 2 is the cross-sectional area 
of the lattice. This is shown as the solid curve in Fig. ||. The fact that the 
data points are significantly above it indicates that the vortex distribution 
indeed behaves as in Eq. (17), supporting our mechanism. 



7. THREE DIMENSIONS 

A superconductor film resembles in many respects the thin three-di- 
mensional model discussed in the previous section, and therefore it might 
seem that the result would apply directly to a real experimental situation. 
However, there is an important difference: In reality the magnetic field lives 
in three dimensions, and the part of it that is outside the superconductor 
also contributes to the free energy. 

For simplicity, we shall now assume that in the symmetric phase, the 
magnetic field behaves as if it was in vacuum. This is not a totally unreason- 
able approximation, because if the superconducting film is thin, the energy 
of a typical field configuration is dominated by the contribution from the 
vacuum outside the film. 

Furthermore, we assume that when the film becomes superconducting, 
it starts to dominate the dynamics of the magnetic field. Essentially this 
means that we can treat the vortices, instead of the magnetic field, as the 
relevant degrees of freedom. 

Again, we have to estimate the typical flux through a circular loop of 
radius £ on the film. The configuration of minimal energy in which the flux 
is <3?, is simply the magnetic dipole, and its energy is roughly 

E($,£)~— (20) 

We can estimate that thermal fluctuations can generate this configuration if 
its energy is less than the temperature, E < kgT. This reasoning leads to 

$ ~ y/nok B Ti = {li Q k B Tif' 2 x $ - (21) 
Dividing this by the area of the loop, we find the number density 
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It is more difficult to estimate £, because it depends on the dynamical 
properties of the superconductor near the transition point, and we don't 
attempt to do that here. Instead, we calculate how large £ must be in a 
realistic experiment so that the correlations between vortices would be clear 
enough to be observed. This requires that the number of vortices inside each 
domain is much greater than one, i.e., 

Ns~^.(»ok B Ti) l/2 >1. (23) 

The critical temperature of, e.g., YBCO superconductors is T = 90 K, and 
this leads to the estimate 

ir 2 h 2 

e » 2 , T « 3 mm, (24) 

which is a very reasonable size from the experimental point of view. 

Experiments like this have not been carried out yet, and there may well 
be technical difficulties involved in the measurement of the spatial distri- 
bution of the vortices after the transition. In particular, one will have to 
be able to freeze the locations of the vortices so that they don't annihilate 
before the measurement has been finished. 

On the other hand, an otherwise very similar experiment has been car- 
ried out recently by Carmi et al.i Instead of measuring the number and 
locations of individual vortices, they measured the total net flux, i.e., the 
number of positive vortices minus the number of negative vortices. The re- 
sult they found was that the flux was consistent with zero, which contradicts 
the predictions of the Kibble mechanism. The mechanism presented in this 
paper cannot explain this result, because it would typically predict an extra 
set of vortices in addition to those formed by the Kibble mechanism. It does 
not make the discrepancy with the experiment any worse either, because our 
mechanism does not actually change the total flux through the film at all, 
and therefore the variation in the total net flux in the final state is simply 
equal to the variation in the initial flux, i.e., very small. 



8. CONCLUSIONS 

I have discussed the scenario of defect formation presented in Ref. |ll| 
from the point of view of superconductors. This scenario is not merely 
a modification to the standard Kibble mechanism, but a totally separate 
mechanism, and in a general phase transition, both of them form vortices. 
Which one of them dominates, depends on the details of the transition. 
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Although estimating the density of vortices formed in the transition by 
either mechanism requires a detailed knowledge about the dynamics of the 
system, it is possible to derive simple and robust predictions for the spatial 
distribution of vortices. Basically, the Kibble mechanism predicts negative 
correlations between vortices, i.e., that nearby vortices have opposite signs, 
whereas the mechanism discussed here predicts positive correlations. These 
predictions can be tested experimentally, if it is possible to measure the 
locations and signs of individual vortices after the phase transition. The 
domain size needed for these vortex-vortex correlations to be significant in 
YBCO superconductors, being of the order of centimeter, is suitable for 
experiments. 
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